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PREFACE 


The author is a Line Officer of the Regular Navy with 
the present rank of Commander and has pursued his graduate 
studies at the U.S. Naval Postgraduate School, Monterey, 
California and Indiana University. I desire to acknowledge 
the sponsorship in these studies of the Chief of Naval Per- 
sonnel, the Chief of Naval Research and the Superintendent, 
U.S. Naval Postgraduate School. 


Y became interested in the theory of non-Newtonian 
fluids during a series of brilitant and inspiring lectures by 
Professor C. Truesdoll at Indiana University in 1951-52 during 
which many of the problems solved herein suggested themselves. 
I served at sea during most of the intervening years and post. 
poned many times the solutions of these problems. 


I wish to thank specifically Professors D. Giibarg, 
formerly of Indiana, now of Stanford University; C. Truesdel!; 
and J.W.T. Youngs, Chairman of the Mathematics Department of 
Indjana University for their encouragement, inspiration and 
patience both at the University and through the mails over the 
years; J.L. Ericksen, Johns Hopkins University, for suggestions 
and review of the draft of this paper; and Dr. J. Weyl, Science 
Director of the Office of Naval Research, who has advised me 
Wisely since my first association with him in 19650. 


I cannot complete this preface without recording acknow- 
ledgement of my beloved wife, Dailey, who has for the past year, 
end from time to time before, coped uncompisiningly with a spouse 
bound to two most exacting tasks, my professional duties on the 
Staff of the Chief of Naval Gperations in Washington and the com- 
pletion of this dissertation, 
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I. OBJECTIVES AND SUMMARY 


This dissertation treats the examination of the effect 
on fluid flow of the theory of non-Newtonian fluids. The gen- 
eral methods of Reiner and Rivlin, as developed by Truesdell, 
are followed. 

Àn investigation is made in Chapter 4. of the func- 
tional dependence of the pressure on the space variables and 
its effect on the possibility of existence of a particular flow. 
This functional dependence is closely related to the form of the 
second order viscosity coefficients. In many cases, making the 
same و و ند‎ of such dependency of pressure as is made in 


the linear case, together with taking the viscosity coefficients 


ee A 
— s 


as constants, forces only the linear solution. This may 


— — 


lead to the erroneous conclusion that the particular flow of 


the non-Newtonian fluid is not possible. This functional 


. dependence is discussed in the general case in Chapter 4. and 


in many of the particular flows studied. 


The stress tensor is developed in such a manner as to 
emphasize the Kelvin and Poynting effects of the non-Newtonian 
theory. 

Several specific flows are considered. These flows 
were chosen az a result of many considerations, the primary one 
of which was, of course, that new results could be obtained. 


Other considerations prompting the selection of specific flows 
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were: (1) the form of the rate of deformation tensor, (2) the 
coordinate system in which the flow could bo considered, (3) 
variation of observable phenomenon, and (4) solvability of 
the equations of motion. 

Chapter 2. reviews the symbols, definitions and equations 
which are further developed in the Appendix. The importance of 
complete understanding of physical —* of tensors and 
their relationships to the other tensor components cannot be 
overemphasized when computations are necessary. 

Chapter 3. reviews the definitions of fluids and develops 
the general form of the extra stress tensor which ís used pri- 
marily herein for the explicit display of the terms contributing 
to the kelvin and Poynting effects. 

Chapter 4. discusses the latter effects and the general 
effect of the pressure on the non-linear theory. 

In Chapter 5., simple Poiseuille flow is studied. By 
means of inversion of series, a general solution is provided. 

An explicit form of the pressure is also produced. It is shown 
that the pressure is dependent on the downstream variable, a 
phenomenon not observed in the linear case. 

Poiseuille flow of a non-Newtonian fluid in a circular 
pipe is completely solved in Chapter 6. In addition, a 


generalization of the Hagen-Poiseuille Efflux Law is given. 
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In Chapter 7., we demonstrate the existence and study 
the structure of shock waves in non-linear continuum theory. The 
steady one-dimensional flow of a trivariate, heat conducting 
fluid with non-linear viscosity is considered in the absence 
of external forces. The discussion is limited to perfect gases. 
Mowever, the proof of existence carries over to fluids with 
more general equations of state. Geometrical properties of 
solutions to systems of differential equations are employed 
to show existence. Comparison is made of thicknessess of shock 
waves in non-linear fluids and linear fluids after choosing 
an arbitrary form of the viscous stress. 

The steady flow of a non-Newtonian fluid in a wedge or 
diverging channel is studied in Chapter 8, This flow provides 
an excellent example of the computational methods necessary in 
the non-linear theory. We are lead to the same elliptic in- 
tegral solution as in the classical case. The pressure dis. 
tribution is, of course, different. 

The problem of Boussinesq is considered in several general 
cases. The effects of pressure dependence are very noticeable in 
this problem. Consideration of the dependency of pressure on the 
distance from the boundary is mandatory. 

The equations of motion of the draining of a vertical 
plate are solved in Chapter 10. Quantitative information is 
provided and a set of curves showing comparisons of thicknesses 
of the linear and non-linear theories with observed results are 


provided. 
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II. SYMBOLS, DEFINITIONS AND EQUATIONS 


General. - Insofar as is possible, standard symbols of modern 
fluid dynamics will be used throughout the text. The symbols 


1 are used wherever general 


and conventions of tensor analysis 
terms or equations are written. Each symbol is defined when 
first introduced. 

In actual computations, physical components are used. 
The beauty and elegance of the invariant form of equations 
unfortunately do not generally persist when the physical 
components are substituted for the flow variables. On the 
other hand, physical interpretation of tensor components is not 
always clear. The equations using physical components and the 
comparable ones employing the invariant forms can indeed be 
quite different in any but a rectangular cartesian coordinate 
system. In fact, it is this difference that lead Truesde11? 
to write "in 1944, 1 first had the puzzling experience of 
getting an apparently wrong form of the equations of elasticity 
in polar coordinates by specializing the perfectly correct 


general tensor equation to the case at hand." 





Ic, E. Weatherburn, Riemannian Geometry and the 
Tensor Calculus. (Cambridge: University Press, 1938). 

2c. Truesdell, "The physical components of vectors and 
tensors,” Zeitschrift ftir ደዎ Mathematik und Mechanik, 
XXXIII, (18953), 345-3506. : | 
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We use tay to mean the physical components of the 
tensor t} 5 

We assume, of course, that a fluid is a continuous 
mediun, with which the reader is presumed to be familiar. 


Equations are derived in the Appendix, but are given below. 


Equation of Continuity. 
Let P be a fluid density 
yt be the contravariant components of velocity 
where ® is the symbol for material differen- 
tiation, 


then the equation of continuity takes the form, 
a 2- 
2£ + (0%), =0 (2-1) 


e 
where "4, " is the symbol for covariant differentiation. 


Cauchy's Laws of Motion. 


Let t? be the stress tensor (symmetric)? 
i 
፥ be an extraneous force field, and, 
x the components of accsleration, 
then, 
6 6 L ۰ ` 2 
በዳ ማፈ ገል | 
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is Cauchy's first law of motion. 
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Energy Equation. 
Let £ be the internal energy density, and 
ዓ be the heat flux vector, 


then 


gi‏ ور حم 


Defining the rate of deformation tensor as usual, 


il 


Y i RU T X. 


t ه‎ » 


A 
and recailing the symmetry of * d we obtain 
: 1 — 
pe = t A ጭ T si (2-3) 
the Fourier-Kirchoff-C. Neumann equation. 


Extra Stress Tensor. 
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ፍቕ 
Defining the extra stress tensor, Vj , by 
GE St (2-4) 
Vis P?it'j 


b 
where 9 is the Kronecker delta, and 中 is any scalar. 


Then the equation of energy can be placed in the form, 
He Veh, 
gon = bizäs + لا‎ ij ~ hi 


| where d is a local statistical parameter (entropy), and 
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For inconpressibie flow dE 9), V 20. In order that we 
may obtain the sane elegant result in the compressible case, we 


choose the Scalar per" and obtain, 
5 -. « 6 | t 
زا 65 م‎ 934 ~F 


a د‎ 92 38 Farm cf the anargy equation. In ረገ flows considered 


——— ——··⸗ 


— ہس TP LS‏ 
.س ا 


© ete ። — — — واو‎ 
وج« دي‎ v ric ~ — — —— —— 141449 —— — —— mn we — — — -— — m ደ መመመ 


፣ መሽ We »)وهه‎ that — stress power is oi RBS Sê, that 8 
ሯረ ٩ 
3 dei; 20 
or equivalently, 
up 
UV. o. 
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6= 9 (en) 
a given forn, 
Fi m KO, (Fourier's Law) 


And the form of the extra stress tensor 
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5 ٣ 2 
where T denotes functional dependence and a,b,c.... are constants 


1 ስዎ” ጨሐቁ፡ — — a 


oí the: fluid, then we have all the equations necessary to study 


5 کس تا 


Pow flow of fluids considered herein and many others. 
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fe. Truesdell, "On thè Viscosity of Fluids According 


to the Kinetic Theory,” Zeitschrift für Physik, CXXXI, (1952), 
273-284, AA c "s 





is used throughout to denote functional dependence 
rather than a specific function. 
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Classification of Fluids and Motions. - 1£ (=0 in all motions, 
then the fluid is said i be incompressible. Furthermore, 11 

e is an absolute constant, the fluid is said to be homogeneous 
and incompressible. 


14 0 « مد ماع‎ all motions, then the fluid is said to 


be piezotropic. 
1f Da pien in all motions, then the fluid is said to 


be trivariate. 


11 € = e (t) , the motion is saíd to be isochoric; if 


0 pit) , isobaric; it P= P(p,t) or e pie t), barotropic. 


All other mptions are said to be baroclinic. 
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111. FLUID DEFINITIONS 


The general form oi the extra stress tensor. - The definition 
of the extra stress tensor is a mathematical formulation of 
Stokes observation! "that the difference between the pressure 
on a plane in a given direction (i.e., a component of stress) 
passing through a point P of a fluid in motion, and the pres- 
sure which would exist in all directions at P (i.e., the 
hydrodynamic pressure), if the fluid were in a state of rel- 
ative equilibrium, depends only on the relative motion of the 
fluid immediately about P, and that the relative motion due 
to motion of rotation may be eliminated without effecting the 
difference of the pressure above mentioned." Thus, following 


2 3,4 5 
Reiner”, Rivlin , and Truesdell , we write 


ti+ ps = Vj - 


where Vi=O, U پا‎ 


ieee =- 
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lo. G. ‘Stokes, "On the theories of the internal fric- 
tion of fluids in motion, and of the equilibrium and motion 
of elastic solids," Transactions of the Cambridge Philosoph- 
ical Society, VIII, (1845), : e ASA 





Zu. Reiner, "A mathematical theory of dílatancy," 
American Journal of Mathematics, LXVII, (1945), 350-362. 


SR. S. Rivlin, "Hydrodynamics of non-Newtonian fluids," 
Nature, CLX, (1947), 611-613. 

An. S. Rivlin, "The hydrodynamics of non-Newtonian 
Fluids I," Proceedings of the Royal Society of London, CXCIII, 
(1948), شک يهو ههو‎ 

ይ 

ር Truesdell, "A new definition of a fluid, I. The 
Stokesian fluid,” Journal de Mathematiques Pures et 
Appliquées, XXIX, (1950), 215-244. 
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We consider only isotropic? non-Maxwellian fluids, and 


set 
١ - زةمه] م‎ + Bc dj + Be, dk Aj] a» 


where ۸ ‚B and T will be defined below. For ease in identi- 
fying Newtonian end non-Newtonian fluids and the higher order 


viscosity effecv:;r, we re-write (3-2) in a more convenient form, 
3 k gi i 
ولا‎ 2۸ 4,34 r2, 4j (3-3) 
v V v Ok 
+ و اه یک + له كا + زوم‎ 


| ۱ 7 
where A+ÉJM is the bulk viscosity , M is the shear viscosity 


and 
Ko = A Go - 4 E 
K A BG, - 24 (3-4) 
IK, -AB'G,. 


(3-3) is an elegant form of the extra stress tensor 


0 


from which we shall take departure in defining all non-Max- 
wellian fluids. From (3-3), the higher order viscosity effects 


become immediately evident. 


Newtonian Pluid. 


If all =© , we have from (3-3) the Newton-Cauchy- 
Y 


Poisson form of the extra stress tensor, 


> e 





Ibid. 223-224. 

a. Truesdell, "The present status of the controversy 
regarding the bulk viscosity of fluids, "Proceedings of the 
Royal Society, CCXXVI, (1954), 59-65. 
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which is the characteristic equation of the classical theory 
of fluids. All results of the theory of linear viscosity, of 


course, stem from this relationship. 


Stokesian Fluid. - Truesde11" characterized the Stokesian 
fluid as a fluid without natural time, i.e., a body which 
offers no response to its previous fora or a "body without 
a memory." | 


If we take 


A= & 
B= 


where P: for incompressible fluids is an unspecified scalar, 





while for compressible fluids, f is the thermodynamic pres- 
sure, and where An is the natural viscosity (dimension = ) ， 
then (3-4) givos 

Ko = Pp Go -ÀA 

Kı = /« 6” ፊ/ 


人 = 2 


which characterize the Stokesian fluid. 


Reiner - Rivlin fluid. ~ We cont 4nuo to follow frseguell. 





m ruesdell, "A new definition of a fluid, 1.," XXIX 
9114. 
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by adding a third material constant Un (dim T ), obtaining & 
fluid with a natural time. The K, characteristics of the 


Reiner-Rivlin fluid are obtained by taking 
እ = 4 
m 
H = tn 


whence, 


Ko = Zen Go -d As 

K, = fen G, سس‎ & fA 

Kr = fu Cn Ga. 

We observe that in the Reiner-~Rivlin fluid, we may 
include Ko in p with no loss of generality. Further, we ob- 
serve that the form of Ki is the same for both the Stokesian 


and Reiner-Rivlin fluids. 


8 وه ما غو شید — — 
- — — ——— — — 
一 一‏ — — — — 


- We consider only isotropic 


—— — sere ee 


The form ot the coefficients G 
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ide and assume Gr TOP TER THén ዳኳ can be placed in the form, 


^e 
= ኋ 3 2 H 
Gr => Gus (BI) (B II) (B III) 
Ll, J, KO 
where 1 ; TI , 11 ; are the elementary symmetric functions 





of the rate of deformation tensor. 

The Geib are dimensionless functions of ,总 
in the case of the Stokesian fluid; and additionally a 
E: in the case of the Reiner-Rivlin fluid. P is simply 


a reference pressure; © a reference temperature. 


We note that 


= e inje , "m! "m - — MILI CANT | 











- "- e bake ir m » , A ‘© ዘ | 
هي‎ E ፡ | Mini, - 6 
22 ጠ ይ 


መ. 

e s oon 

Abe dew e A 

deem e A 

sw ey «‹8/ዘ!ነ ሠ/(1/ዕ››»፡፡‹ መ 4 
7 eee, د مت مه لها‎ lg Er ې‎ 


ey ana سم عم‎ par nu www Le 
dem i Jm 


abe የመ ww — E LS 
ne rcrum vium 


۱ 


* 2i BID A) (za) ve E, 5 el 





13 
A= AB acı 
(3-5) 
سا‎ ۸6۰٠۰ 


in both the Reiner-Rivlin.and Stokesian fluids. 


er 
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IV. GENERAL EFFECTS OF NON-LINEAR THEORY 


The extra stress tensor in extensio. - In order to examine 
clearly the general effects of the non-linear theory, we vrite 
the components of the extra stress tensor in extensio, thus, 


from (323) ; 


vi = Ad + 2pd, +K, +K + Kı dud, 
Yi بت‎ 4 ሷኒ + ره‎ 十 Ke 二 KR + Ki AR 


Vi = A A+ 2 P- d, + Ko +K, As F Ka Ah à (4-1) 


AE 2p + +K Ai 
Va == «p A, +K, &, +K, A, a 
a, = ዴዶ Az +K, A +K 0% 


Kelvin Effect. - The terms Ko tK, d, Kot KG, ana Ks t ۵ 
of (4-1) represent the Kelvin effect, a phenomenon that occurs 
in the theory of non-Newtonian fluids GA is not ችንን in 
the classical theory. If a hydrostatic pressure equal to the 
negative of Ko is wanting, the fluid will tend to change its 
volume. In tbe case of an incompressible fluid, an arbitrary 
hydrostatic pressure may be added to the stress and no new ef- 
fect is noted over the classical case. This effect has been 
termed the Kelvin effect by Truesdell} after Kelvin who pre- 


dicted it theoretically”. 


Libia. 


Ze. Thomson (Kelvin) and P.G. Tait, Treatise on 


Natural ሩፌ ተፅ (2nd ed. 2nd half; Cambridge: University 
1 و655‎ 
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Poynting Effect. - Each of the components V, ሄ Tu and vi 

of the extra stress tensor contains a term of the form K: 2,4; 
etc. Hence normal tensions must be supplied in the planes of 
motion; if such tensions are wanting, the material will tend 
to change its length. This phenomenon occurs in the theory of 
non-Newtonian fluids, without regard for compressibility and 
we call this the Poynting effect, following Truesdell?. 
Poynting first noted this N =. Ha stated that "when 
a wire is sufficiently loaded to be straightened, it is leng- 
thened by twisting by an amount proportional to the square oí 
the twist and with a given number of turns, inversely as the 


` 5 
length." Poynting later reported the results of measurements 


of diameter contraction, which verified his predictions. 


The effect of pressure on the non-linear theory. - Let us 


write the extra tensor tensor in the forn, 
l V ? x 2€ 24 
- ود دا‎ th Aj +E Ma! 


where the Fe are defined by the corresponding coefficients 


in (3-2); then the general equations of motion of the steady 
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rruesdoll, "A new definition of a fluid. 1.," XXIX. 
4 
J. H. Poynting, "On pressure perpindicular to the 
shear planes in finite pure stress and the lengthening of 
loaded wires when twisted." Proceedings of the Royal Soctety, 
LLXXII, (1909), 546.599. | 





SNR. Poynting, "On the changes in the dimensions of 
a steel wire when twisted and the pressure of distortional 
waves in steel." Proceedings of the Royal Society, LXXXVI, 
(1912), 534-561. 
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16 
state flow of an incompressible fluid may be written in the 


form 


u ` — E 
0 ካረ ل‎ ጨጨ 9; + = ፉ | ۷ 
et; SE SM -- à f, à; A, Xi, A; — (4-2) 
Ap 
Let 2 be constant and lot «4; and ۶ be the velocity 6 
pressure which satisfy the corresponding Navier-Stokes equation, 


that is, let 
"D 2n 11 
eh دم - ]نع‎ + 7 ወጀ 
We briefly discuss the conditions under which W; also satisfies 
(4-2). To this end we divide the pressure in the latter case 


intc two parts, thus Pf: P . Then 4, will satisfy (4-2) 


whenever 


Pit D (4t 4t], , . (4-3) 


Generally, the right side of (4-3) may or may not be a gradient. 
When it is, we can obtain solutions of this type. 

Three possibilities occur from (4-3): 1) 7, =O , in which 
case P: f and the flow is Newtonian; 2) there exists a non- 
trivial and additional pressure gradient, P, y ‚ which satisfies 
(4-3) and must be added to P in order for uU, to be a solution 
also in the non-linear case; and 3) there does not exist a 
function, ۳ , Satisfying (4-3). In this case, secondary flows 
may have to be introduced in order to solve the non-linear flow, 
if, indeed, solutions exist at all in a particular case. 


Yt is further noted that a special case of 2) above occurs 


when (ፈኒ A: ),,*0 : 
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V. SIMPLE POISEUILLE FLOW 


General. - This flow was first introduced by Euler in 1757, 
It is considered herein to demonstrate the use of inversion 
of series to obtain solutions. 

Let there be two infinite parallel planes, the lower 
one fixed, the upper moving parallel to itself at a steady 
Speed, and assume that the fluid adheres to the boundary. We 
take the X «axis in the direction of motion of the upper plane 
and the a Axis perpindicular to the planes. Steady state is 
assumed. 

Then 

XG) = X= (909, 0,0) xG)20 
e iw e| 

d (i$) = iq 9 o 
0 ጋ © 


with elementary symmetric functions 


Te De, Dakar 


and 


۰( = 31خ -) + 2 میم 


The physical components of the stress tensor are given by 
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The Kelvin and Poynting effects are observed even in the 


— — — جو ںہ کا س سے لے ن سے ود مسد > 


case of isochoricflow. A pressure must be supplied perpin- 
dicular to the planes; if this pressure be wanting, the upper 


plane will rise or fall, depending on the sign of h : 


The equations of motion. - The equations of motion are 
D e 3 (- p+ Ko + £3 &)*- $; G gez) 
=- 2 o ። (5-1 
0 yv q [K. + 244)) gës و )کو‎ Ka) 
ጋ 
0= 22. 
The solution. - From (5-1) 2, we have 


296 — o 
In du l 


hence, from )5-1( 


2 $) ME. (5-2) 
y 


where, 


F=ABG,=K,+2pm. 
Integrating (5-2), we obtain 


Ü ua = 49% = Ày +A, 99 


the resistance per area of the flow. 


' 
Again ያ፦ኙ05 (5-31) 1, : 


= AX + P (y) 


x 
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19 
and from (5-1), 


-p t^ La د‎ = -Å t Pin). 
Hence 


$ lu) = (zz) = A = ዚር )-+ ره و‎ 9789 


wherein Ay are constants of integration. 
Rivlin congidered the compressible case and, further 


assumed A, =O , obtaining the resistance per area 


፡ሪ/።ፊ -#ቆ ሩጋ 


ፒ ረ/ዶ) = @, g' + 0,9" + Ay 4""ፉ - Qc s. 


has pointed out that, depending on bow many‏ ا 
terms we take in this expansion, there will be an odd (or‏ 
infinite) number of functions + (i.e., a certain velocity‏ 
profile) which will yield the same resistance per area.‏ 
However, we shall assume that 3 is an analytic function and‏ 
furthermore that Eliz) is a monotonic function of g’. This‏ 
will permit inversion oí the above series, obtaining 9' us‏ 
a function of 4 Ye then integrate term by term to obtaín‏ 
the function 34) . This will be accomplished shortly.‏ 


The rectilinear shearing flow of a compressible 


1 : m nn nn 
Rivlin, "The hydrodynamics ......," CXCIII. 


D E "A new definition of a fluid, I.," XXIX 
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non-Newtonian fluid was first suggested by Reiner” and 
Rivlin* and elaborated upon by Truesdell". Further work by 
Trussdeli? treats the insufficiency of the shearing forces 

in the non-linear case. We also considered the rectilinear 
Shearing flow of a compressible Stokesian fluid and reported 
"the possible velocity profiles ና now depend upon the pressure 
P as well as the resistance" and that "as the pressure is 
reduced, the resistance is diminished as should be expected 
in a gas. Eoth in the classical theory and in the second 
approximation, we have ۸۰5 وم م7‎ * Ay, but in ማት third approx- 
imation, tho resistance is diminished by سو و‎ 

the value predicted by the classical theory for the same 
velocity profile Y" 

Truesdell's statements regarding the possible profiles 
of 7 are certainly correct. However, since we are dealing 
with physical phenomena, one would expect. a more unique be- 
havior of the function a 


Returning to (5-3), we rewrite, 
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Be Truesdell, "The mechanical foundations of elas- 
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4 7 = -Ay th. (5-3) 


where the left hand side is a power series in y. From the 

assumptions regarding analyticity and monotonicity, it is 

possible to invert (5-3) to obtain 3'as a power series in Y. 
In either the Reiner-Rivlin or tbe Stokesian fluids, 


we have 


1, = A BG, = MaG, 


where 


> e J 
/ 2 
(3, => 6 ۳ Ly) | . 
了 2 0 
We assume the G4 ojo constants. 


Then 
p ao, B A Bas pt 44 p 
1 ጠጋ fen (2. መጫ E d + à 2 9 — 9 ge e 2 7 


and from (5-3) and (3-5), 


Ma a Bs‏ موک رو ونا 
AAA A?‏ 
A pee (5-5)‏ = 
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22 
Reverting 5 (8-5), we obtain 
9 (y = — zr (A y- A) 一 = (ዳ.3- A.) 


رد ٢م‏ 0 Ž y‏ ر بے 
ራፉ. (My -AN‏ 


Whereupon, integrating, term by term, 





Aly) = Asty ይዮ፦ La (gan jo >> E 


Again we observe the effects of assumptions in the 
classical case. For example, if A, is assumed to be zero, 
we immediately obtain a linear profile. 

In (5-6), the condition of adherence to the bottom 
plene gives Á,*O; that of adherence to the top plane will 
yield a value of ۸ or A x. A further boundary condition is 
necessary, Bay, g” at the upper plane, in order to determine 


the profile uniquely. 


The pressure. - The pressure was determined in (5-4), namely, 


# = Ko +9 7% - ۱۷ 


In the case of an incompressible fluid, Ke can be dis- 
regarded as previously discussed. In the case of a compress- 


ible fluid, Ke determines the Kelvin effect. In both cases, it . 
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observed that the pressure is linearly dependent on the 
downstream variable, 
| Finally, in tho case of isochoric rectilinear sheer of 
Sa COmpressible Stckesian fluid, terms of the form B*"= (c^ 
which appear in the explicit expression for q (4) in (5-6), 
indicate, as Truesdell has often stressed, the appropriateness 
of the Stokesian theory for low pressure phenomena in gases. 
in this case, there is no longer a general explicit expression 
for the pressure since in general, ይ s A and 3 are power 
serles in de . However, in the case that fo and m are 
constant, 
۸ Ak E 
Tu sd 


and the expression for the pressure remains valid. 
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vi. POISEUILLE FLOW IN A CIRCULAR PIPE 


General, - We briefly consider laminar flow through a cylindrical 
Dipe under steady state conditions in order to derive an extension 
of the classical efflux law. We take the axis of the pipe as the 
Z -axis in a cylindrical coordinate system 2,6, 2 . Adherence 


to the boundary is assumed. The velocity field is giver by 
X(t) = ረ2. 9. #“/ 
and the acceleration is zero. 


The equations of motion. - The equations of motion are given 


by 2 La رو‎ 2 

on Sere + eh lenrg 4) 

. 2t. 
5ያ , 3 (tq Hf). 6- 

ez Ser Selz? 7 T 
The form of the nressure. - From (6-1)5, we deduce 

5ዎ ta) 

Oz. 
or 


p= 2 ዜ (ዱ) + ፳/) 


and therefore 


iP - ፖ ክ(ጾ) + ፳ ር) 
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but from (6-1), we note that 


2t = pa) 
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/ 
hence kh =o, h = A, , and the pressure gradient is there- 
fore constant downstream. 


Again from (6-1)5, 
+ 9۳: +چ-)‎ FP EAR )6-2( 


Integrating (6-1),, by parts, we obtain 


: 3 
د م2 - مر‎ -Azr sy hr a da 
and substituting for g^ from (6-2), we obtain 
r 
P = Fe ~ zg age, ፦ 45 ja E da. (6-3) 


The flux, - The main flux through the tube of an incompressible 


fluid is given by 
مت‎ 
03 1 2 7۳ 2 
٥ 
Integrating by parts, we obtain 
la 
2 a ng 
IL 
Q. = ame | (541, | ERA 
The first term in the bracket vanishes because of the condition 


of adherence, and therefore 


* rp 
a - ame, | ጅሙ 
If ዳ is constant, then from (6-2), we may write 
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a parabola, the classical resuit of Poiseuille. 


Further, from (6-3) 


SAh 
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Defining the inlet pressure, 


زو 8:8 


and the outlet presgure 


bip 


where 1 is the length of the pipe, 








e <= -A Lt: 


we have 


# ኤሪ 
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top 
A, = ki 
and, tne flux 
7 | (6.4) 


This expression for the flux was first given experimentally 


"4 
by Poiseuiile Pers and is termed the Hegen-Poiseuille Law of 
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J.1..K. Poiseniile, "Experimental research on the motion 
in tubes of very mali oinueter," Academie des Sciences Comptes 
Rendus, XI, (1840), 9801-967, 1041-1048, 
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“Ibid, XII, (1841), 112-115. 
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Efflux from a pipe. The formal derivation is due to van 
A generalization to this law is derived below for the non- 
linear case. 

We observe that the velocity distribution and the 
mass flux are unchanged from the classical case in the second 


order ( ورگ‎ constant) theory. 


The solution by inversion. - We return to the dilfereniiul 


— — — 
nn ——‏ دیص همد 


equation, (55), ® 
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፡8'8:(፡8) = - Ar. (6-2) 


By methods similar to those in the previous chapter, we 


. derive 


3 pa 
fin Ai ساب‎ (as 


ቋ‹) ብቃ zu (ee! + 32 


The mass flux in the case of a non-Newtonlar incom- 


pressible fluid is given by 


a= — 
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س 
Computing only tbe first two terms, and taking‏ 
A= EXA , ss before, we obtain‏ 
é‏ 
ያ —‏ 5 


an extension of the Hagen-Poiseuille Efflux Law, which is 
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more conveniently written 5 | 

6.= ፪፡. [(ዶ(ይ-ሥ)) ፉ for Bia. (see 

= * AE ራሥ ی‎ | . 

Poiseuille used a "pipe" of radius ranging from 0.007 
to 0,325 mm. The new law predicts that if similar experi. 
ments were performed with NC sufficiently large, yet 
the flow still laminar, that a measure of the efflux would be 
substantially different from that given in the classical 
efflux law. 

In general, gr) will be an even power series in 
Ers. and Q will be an odd power series. The classical 


£ 
Hagen-Poiseuille Efflux Law is valid whenever 


pA. died. 
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VII. STRUCTURE OF SHOCK WAVES 


General. - The existence and structure of shock waves in the 
non-linear continuum theory was established for the fluids 
described herein jointly by Gilbarg and the author”. In 
addition therein, we computed shock thicknesses in monotonic 
gases and compared these results with those of Thomas“ (con- 
tinuum theory) and Mott-Smith" and Zoller* (kinetic theory); 
and further demonstrated that the besic methods of extending 
the theory to polyatomic gases remains unchanged in the con. 
tinuum theory provided that fluid moduli can be determined. 

That portion of this joint work pertaining to the struc- 
ture of shock waves in the non-linear theory is included as it 
relates to the topic of this dissertation. 

The steady one-dimensional flow of a trivariate, heat 
conducting fluid with non-linear viscosity is considered in the 


absence of external forces. The flow approaches finite limit 


values at X=3*00 . Such fiows display the character of a shock 








lp. Gilbarg & D. Paolucci, "The structure of shock waves 


in the continuum theory cf fluids," Journal of Rational Mechanics 
and Analysis, II, (1953), 617-642). 

- Thomas, "Note on Becker's theory of the shockfront," 
Journal of Chemical Physics, XII. (1944), 449.452. 
"Wn. M. Mott-Smith, "The solution of the Boltzmann equation 
for a shock wave," Physical Review, LXXXII, (1951), 256-274, 


Ze. Zoller, "On the structure of shock waves", Zeitshrift 
für Physik, CXXX, (1951), 1-38. 
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wave in that they differ from their end states at X = عد‎ o0 
only in & small interval of rapid transition. 


The velocity field is given by X su and the rate of 


deformation as $+ Tu ; Where the prime will denote differ- 


rm مود‎ o — — — 


entiation with respect to X throughout this chapter. 


mat up — — — 
u نموت‎ ው 4ም= mamie x 
— — 
— — — 


Equations of motion. - The equation of continuity, Cauchy's 
first law, the energy equation, Fourier's Law and the extra 
stress, respectively, take the form 
/ 
(eu) = 
puw=t 
? / : 
00.6 =Fu’-g 
^ 
: -X 9 
y= Cr 


Elimination of f and $ from (7-1) and integration of the 


(7-1) 


resultant system yields 
gut Me 
gurp-vV=P 
eu[£rerz]-vu-Xxo'=eE 


(7-2) 


where My, P and E are constants of integration (7-2) ex- 
presses, respectively, the conservation of mass, momentum and 


energy. Eliminating ٣ and redefining 
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2 ርቀ Sen C= Am? 
we obtain 
p +blu-a) 
pat ሪ WS N (7-3) 
ير‎ ©“ > h fE- أن - *زيه-ه) عل‎ 


The form of the viscous stress. = We propose to study the ef. 
fect of uon-Linear viscosity on the theory of the shock profile, 
an application naturally suggested by the large velocity grad- 
ient present in tne shock front. 8 
Of the possible generalizations to fluids obeying a non- 
linear stress law, we select one in which the viscous stress 
may be a non-linear function of the velocity gradient, but 
remains independent of temperature gradients. To fix our ideas, 
we choose the specific volume, , and the temperature, © , 
as variables of state and define the viscous stress ሸ፲ =ህ (u',v,0) 
with the following properties: | 
a. Y(0,v,0) =0 for all ~, © 
b. 29 (6) 5) = ሥረጊ/፡ ይ < 0 rer ail A, و‎ 
Property a. asserts the essential characteristic of fluids that 
the viscous stress vanishes when there is no velocity gradient. 
Property b. states that for small velocity gradients, the stress 
15 approximately given by the linear law and the fluid essentially 
becomes ge classicui one. 
| For convenience, we consider the function y = y lu, e, 6) 
which we identify with V (u', v, 8) by means of (7- =2) 1 and re- 


write (7.3) 
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4 6“ - t [t- t0-2-c6]. 


(7.4) 


For convenience in exposition, but without essential 
loss of generality, we limit our discussion to perfect gases. 


However, the same proof carries over as weli to fluids cbeying 


the general equations of state considered by Gilbarg” in 


proofinr existence in the linear case. For these gases, 
p=Ro£, E=CyG6, R= Ge tr 


where Cy , Cy, are the specific heats per unit mass at constant 


pressure and volume, respectively. We set r-z , as usual, 


then define 9 > رس ول‎ 


6 


Following Becker”, we introduce dimensionless variables 


by means of the substitutions 


I w> RO 
د‎ /ፆ-ቴ- 6-፪ጀሸሥ 


and also define 


E m” = 
ሽራ - ; OU - XE 


and the function 





9p. Giibarg, "The existence and limit behavior of the 
one-dimensional shock layer," American Journal of Mathematics, 
LXXII, (1951), 256-274. 
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中 = VU (ሠረሠ,ሠ ፀጋ, ulw, w, E), O (w, w, eil 
= ሀየ ረሠ" ሠ, ወ3. 
After the &bove substitutions, (7-4) becomes 
pl, wW, Ô} = رن‎ + 2 -/= M(w,0) 


x € = 6 - 3 |) - “(س‎ + <1 = | (ሠ. 5) E 


6 
1 
1 


0 carries over the properties oí V (o e, 2) namely 


a'. له‎ (oe, وله‎ 5) = o for all «, Ó 
bit. 22 (0,w,8)* A (w,ö) 20 for all oi, a. 
0 


By the implicit function theorem, these properties have as 
an immediate consequence the existence of a function Y (a, O) 


such that 
W= w'(M(w4,5), W, 5) = y (w, 5) 


in *&ome neighborhood of the paraboia M-Ó, and such that 
¥ =0 on this curve. Naturally tbe extent of the domain of 
definition of Y (w, 2) cannot be deduced from the purely local 
properties such as a' and b', 

In the W-@ plane (Z-plane), iet the two parabolas, 
M0 and L=20 intersect in the points ፲ =(..8), 7፻(ሠ.,, 6) 
with وله‎ < 0J, . We see that such points could correspond to 
initial and final states of a shock wave, since Zo and و‎ 


would simultaneously satisfy (7-2) and therefore satisfy the 


Rankine-Hugoniot shock conditions, namely, 
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where the values of kao, &, , etc are obtained in the obvious 


manner. 
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Equivalent conditions for existence. ~ We will now reduce 
(7-5) to a form —— emus techniques developed by Gilbarg/ 
may be applied. | 
| We make the further assumption that y (4,6) is defined 
throughout the region, R, (see figure 1.) bounded by the 
parabolas. 

Since #= 0 implies W’= ©, which in turn implies ^f-O 
(property a'), we observe that ¥% =O only on the parabola 
M ረጩ, ቋያ =O . Therefore Ylw, 6) is of one sign throughout R. 


Furthermore this sign is negative since 
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joining the singular point Ze and Z,, then clearly we have 


| succeeded in proving it for the system (7-5). 


Proof of existence. - The proof, patterned after the one for 
the linear case p. rests on tne following facts: 

(i) Zo and Z, are, respectively, node and saddle point 
of (7-6); 

(ii) there 15 5 unique integral curve S(x) of (7.6) 
that approaches the saddle, Z, , from within R as X> œ, 

(iii) © has a negative slope in R and cannot intersect 
the boundary of R. 

(iv) S cannot terminste in the interior of R and there- 
fore must approach Zo. 

Steps (iii) and (iv) follow, as in the linear case, from 
our knowledge of the signs of M(w,8) and L(w,Ó) in R. 
The parabolas divide the neighborhood of Z, into four regions 
0455 which the slopes of the direction field are alternately. 
positive and negative, and (see figure 1.) in R the slope is 
negative, and in the direction of the second quadrant since 
Mw, 5) «oO S, L(w,5)>0 throughout R. We observe now by con- 
sidering the signs of L(4,9) ana M(w,8) that on the boundary of 
R, all integrai curves are directed out of the region as X in- 
creases and therefore into the region as x decreases. Hence if we 


follow the integral curve S(x) backwards from Z, in the direction 


vs - 
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of decreasing X (and increasing له‎ ), we see that it can never 

leave R and must approach the other singular point sS ጂመቅ -ወዐጋ. 
It therefore remains only to discuss the nature of the 

Singular points at Z, and Z, and thereby to establish (i) and 


(ii). The characteristic equation of (7-6) at Z, and Z, is 
| م۳‎ - © Ys 
=) p DA e? 
| 6 aJ او‎ ረ 2۸2 502 


We observe that 
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given in the linear case, that is, 
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By inserting the values of M and L and subsequently 


evaluating at /, and Z, , and after simple but tedious 
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algebraic computation, it is not difficult to see that the 
roots of the characteristic equation are real and cf opposite 
sign at Z, and real and of the same sign at Lo. Thus Z, is a 
saddle point and Z, and unstable node. There are exactly two 
integral curves that approach the saddle 2,25 x-yco0 and exactly 
two that approach as X>-&0 . These pairs correspond to the 
negative and positive roots, respectively of the characteristic 
equation. Two members of each pair have the same slope, but 
approach fron opposite direction. Thus (1) and (ii) have been 


demonstrated and the existence established. 


Numerical computation, ~ The numerical calculation of shock 
thicknesses in the linear fluid under discussion requires 
somewhat exact knowledge of the viscous stress, y (44,6) . In 
an effort to arrive at an estimate of the infiuence of non- 
linear viscosity on the shock profile, we consider a gas for 


which we assume a quadratic viscous stress of the form 
/ p 
: Y = — au په بپ‎ -٨ر‎ 6e) uw )7-7( 


DPE At is the ususl shear viscosity and A (v,0) 45 a co- 
efficient still to be determined. For lack of specific empirical 
information concerning A , we arbitrarily take a value suggested 
by the Burnnett approximation of the kinetic theory. In this 
approximation, the expression for viscous stress in one-dimensional 
form in the absence of external forces and which contribute to (7.7) 


is S 
4 , 2 /4- (b 
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where BD, , D, , Wy, are dimensionless quantities whose values 
depend on the molecular model. 
Simply to serve as an example, we take the values for 


D,, 0, and D, given by Chapman and Cowling!! 


4 
دوس مون لد وتوو نا 


Substituting our arbitrary expressions for Y into system (7-5) 


|+ + 一 | 


Z- Áu 
e 


Oo = 8 -5 ርነ-ሠነ” +ሟ] 


The positive square root is taken for w’ in order that W'=o 


we obtain 





je lj’ = 


imply M=0. For purpose of calculating the shock profile, this 
System can now be treated as was that resulting from the Navier- 
Stokes equation, If ww ei? ; then the expression for 4 is 
simply م‎ 45 T : 
In order to compare with previous results obtained from 
the linear theory, we specialize to the case of Helium for which 
= 647 and eL 11 we designate by PAY and Aguad. the 
shock thickness according to the linear and non-iinear theory, 
respectively, our sampie calculation show: 
at Mo =2 qud. 6 
him. 
guad. _ 130 
and at ا‎ SCH : 


«ሪኡ .مدو(‎ 
a 


"t Chapman & T. Cowling, The Mathematical Theory of 


Non-uniform Gases. (Cambridge: University Press, 1939) 
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Thus the shock thickness increases in going from the linear 
to the non-linear, the percentage increase being greater at 


the higher Mach numbers. Of course, the significance of these 


numerical results is relative, since the choice of the form of 


the viscous stress was arbitrary. In a specific real fluid the 
effect of non-linearity might be greater or smaller depending 
on the empirical value of the coefficient of the quadratic 
term. However, the discussion shows that the potentialities 

of the continuum theory are by no means limited to the Navier- 


Stokes Zquations. 
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VIII. STEADY FLOW BETWEEN NON-PARALLEL WALLS 


Introduction. - We consider the laminar flow of an incom- 
pressible fluid between non-parallel walls, under steady 
state conditions. Se take the intersection of the walls as 
the Z -axis in a cylindrical polar coordinate system and as- 
sume adherence to the walls. The coordinate, 7t, , the distance 
from the Z-axis is assumed greater than zero. The equations 


of the walls are therefore given by © طء‎ # 


General. - Radial flow and absence of secondary flow is assumed. 
x) = (u,o,o). 

Then the equation of continuity takes the form 
د‎ መ. =0 

which ig satisfied by 
۲۲ 4» 2 9, )9,2( 


anc, finally, if we assume that * zO, then 


ዜፌ = 8:9). 
ما أ‎ 
The inertial terms are therefore given by 


X(i)= (-%-.0,0) 


y 


and the rate of deformation tensor as 
-2% ዓ 
Aii =- | ود و‎ 0 
5 = 


AY 
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with elementary symmetric functions 


18 =o ]إل‎ = - legg). 


The stress tensor is written 


ም -- 1 © 6 M | -29 y o 
1 = سه‎ 0 — o 
qe > P > ale E | ሆሪ 


K, [2% Y O |م‎ ٠ © .9 
سس‎ | 8۵ 270 ት ታታ 0 43+3 O 
۸۲۰۱ 6۵ oe ۱ 0 0 0 
The Poynting effect is observed. As in the case of flow in 
a pipe, Tür) differs from the pressure by a quantity which 
is always greater on ihe edges. Hence, if this pressure is 


wanting and the walls are terminated at ቪሺ= hb ,رح«‎ then the 


emerging fluid wili swell. 


Equations of motion. 


Let 6 «IL. 


Then 


ds d (97 


and the stress components may be written 
E(u) = و ټون ې‎ + 0 
b(a2)= -p ++ 70 


#9 
CUR Qu 


C (33) 


1 


0 





7 (3) = L[ (23) = 0- 


. a 1 
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The equations of motion are 


ı 2(9ኾ) 2(#6) ዱ:2 20 


eg مهد‎ 
ፕሮ: Bree wore In ፳፻ وه‎ 
Y 1 ኃ( F, 206 7و)2 سل‎ = 
o = اطع د لقو ج + غو‎ en 


0 <= ።: መመ se . 
2 
Differentiating (8-1), with respect to ቀ and (8-1) with 


respect to FL, and subtracting, we obtain 
E ም — 2 (v1 et 3 (v 7) 
4937 = —4n ra + 5 ጣቸ 
5 3(ጳ.ቾ MIE) | 
mido + n gi 


In order to separate explicitly the dependence of the 


flow variables on the space variables, we introduce several 


substitutions. We first take 


_ Hd) <- 2 " " 
6 = 4. n* 2 Hid) = 44 Th (8-3) 
Then , 
290 . 46 ae . 200 
orn " FE or ۴ 
58 _ dp 4 ሇ / (8-4) 
29 = 9(-”ዛፒፒ. -Sp > 0H] (Ln 41) | 


ፓወ ዱ ; 
26.29 ጨ 一 x 9 (lu 541). 


Substituting (8-3) and (8-4) in (8-2) and after some manip- 





ulation and gathering of terms, we obtain 


مه 


pr 


zg 


AA 


4 ووم‎ - Z'^e'[s (i. H)" — ieg + 169 (LH) 
tE OJAS H + gm HH) 209 (La) 83] 


)8-8( و + و۱4 + 


independent of f. : 
Substituting for © 


6 ه -ووم + + ره و + هه ዕኔ እ 8(9) + 255 እድ‏ 


and rewriting, we have 


where 


B5) =[ YH” + 163 HH —163'H) 
P(9) =[-23"H-—3'H"—203H'+ 89'H] 
Polo) =[ 3" +۷ 


af 
~ Suppose now that we take +, 
وی‎ 
» 
= د‎ Ay O 


»20 
are constant. Then 


Solution. of the forn 


where the Q,s 


a < = 
A9 2. d 


di, -ቃ ሃው ره‎ ወ” ځ‎ 


B... that ^ S ¿e and substituting for Fy and 


its derivatives in (8-6), we obtain 
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or equivalently, 


O = > poar -»PH+ RH Ja» —— 


(8-7) 


ተ(8 -5 ዘ) ጁ% ተ 5 5» + ቶተደ2# 


Equating the coefficients of the powers of p on each 


side of (8-7) we obtain the following equations, 
Pe@. + 4043 =0 (8-8) 
[:(9-0 &,- » P.H + Pe Hî] as = 


۱ fowm Xal, 2,3, enet (8-9) 
Supposing, in (8-9), that Qy%O for all ۷2 | 
then 
Y(%-1) BR - YP. H + Po HT = O for al1 ۰۱۷۸ | (8-10) 
is a system of linear algebraic equations in the three 
unknownsp. » P,ti> Pot 
For any (X, P, y , all unequal, the determinant of 


coefficients of the system (8-10) is given by 


& (471) - ۱ ol 一 以 ) 
6(8-0 -5 | = |ø -Æ | 
YLr-/) -Y i yY? -yY و‎ 


the Vandermonde determinant, which has the value 


(A-B) 8-1) (7-4) 


and vanishes if and only if at least two of 4 | 2 ; ጋግ are 
are equal, contrary to hypothesis. Hence it cannot be the 


case that all (Ly are non-zero. 
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Repetition of the above argument demonstrates that all Gy?! 
are zero except lor at most two different, say and 6 ፡ 
In this case (8-10) reduces to 
Aqai) Pa -APH =- Pa H> (8-11) 
م‎ )» -١( ثم - و6‎ PR = - PH" 


or 


_ Pe” 
Pa = a 


P, = —— = ዴ/ያ/ፆ።ፉ“«-፦// 


where )ني‎ and ዶ are to be determined. 


Now from (8-8), since ($52 2 ir 


(8-12) 


q” +49" + E 22 =። 


which is the differential equation of the classical case 


Y 


l 2 3 4 
Stuidied by Jeffrey , Hamel, Harrison and Oseen , and which 


after integration, ylelds 


8” - ln - ur 一 Zi SCH (8-13) 


where A and ደ; are constants of integration. This equation is 


E E 


Je 7 Jeffrey, "The two-dimensional steady motion of 


a viscous fiuid," Philosophical Magazine, XXIX, (1915), 455-465. 








e 
"G. Hamel, "Special motion of viscous fluid," Deutsche 
Mathematiker Vereinigung, XXY, (1917), 35-60. 


a Harrison, "The pressure in a viscous liquid moving 
through a channel with diverging boundaries," Proceedings of 
of the arme Society, XIX, (1920), 307-312. 


3 ር W. Oseen, "Exact solutions of hydrodynamic dif- 
ferential equations," Arkiv for Matematik, Stockholm, No. 24, 
(01927). No. 22, (1928). ۱ 
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further studied by nenne of elliptical integrals. 
Ve write, in lieu of (8-15), 


| * 
3 = — [2E Gear mice gn 


where two of the G,, Ge , Gs , are constants of integration 


and 
د‎ 8 
C, مع‎ + Ey dies = 
e 
kis requires at le&st one of the Êr should have a negative 


resi part. XZ we write 


Ra ce. ) ቃ Qu. (6x) » Qu Ces} 


then 
Ga (Ey) € — Ze 


the sign of equality being valid only when the Cy all have 
-¢he sume "1 part. The different possible cases have been 
discussed by Hanel. 

Returning to (8-11), if we substitute the value of 
from (8-8), we huve two equations in 2: its derivatives and 
their powers.  (S-11) is then a sot of diaphantine equations 
in Á end P . 1 have not been abie to determine whether non- 
zero solutions exist in this caso. ፲፻ such solutions do not 
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collecting terms, we obtain 





2 (p-ho) 2 L 35 (p-20) = ep 


0 FL 
where 
Q, (>) SMG + hpg' tle (8-14) 
We may therefore assume that 
2 * = 
r 3 (ዞ፦”ዓ = ደ ወ,(ቃ) (8--15) 
and 
2 ` 5 A (8-15) 
r$, የ 1.9) = Qs (9) 
with 
Integrating (8-15)7 9 with respect to ft and ቀ , we 
obtain 


Y + ke, (¢)‏ - = 6 و مد 


p-e = $ (asah thal 


since SÉ o from (3-1)3. 


(8-17) 


Integrating (S-12), we obtain the second derivative 


of the function qd: 


I^ 9" 9 43093 + peh, = O 

which, when combined with (3-14) and (8-15) yields 
ደ 9, (ቀ) + (33 (ቀ) - قا( ء‎ d Zä ë + MA = O- 
Supposing now that we take, for convenience, 
QLI — 43 — ae 


— / 
and Q3 = R p 
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49 
then (8-17) become 


p -天 = = + pe + kı($) 
ች - - 4,0 


and finally, by comparison 


" مود‎ 360) ፉ ሥብ. — O + consi. 
f m an 


ZEE + kala) 


The above analysis demonstrates that the flow of a 
non-Newtonian fluid in a wedge yields a velocity profile 
identical to the flow of a Newtonian fluid provided no 
secondary flow exists. The pressure, however differs frozi 
the classical case by an amount 7 O . It is noted that 6 
restriction has been placed on ዴ , while F, was assumed Cen- 


stant in arriving at this conclusion, 
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IX. THE PROBLEM OF BOUSSINESQ 


General. - in 1830, Boussinesq! posed and discussed the 
problem of "how tne retarding influence of viscosity trans- 
fers itself throughout the entire fluid mass, beginning with 
the boundary, as soon as the state of rest ceases." 
Boussinesq postulated that "at a given time, a constant force, 
. parallel to the horizontal exis“ is exerted on a unit mass 
throughout the fluid, in a manner to set the fluid in motion, 


except that the bottom maintains its own adberence to the 





boundary." 


* 


We therefore take x” as a set of rectangular cartesian 
coordinates and consider an infinite fiuid occupying all space 
in the upper half plane, X23 O , and take the Z-axis vertical. 


Then, for 

t< O; se 02 و‎ ጁሄ*= ዕ 

ደዕፕ 

t20; «۲ )۵( =(k,0,0) 
X*=Xx() = (u, 0,0) 
w = g(t). 


The boundary conditions demand that 
a(o,t) = © 7 3 (2,0) 7 o. 


ly. 3. Boussinesg, "On the manner in which viscosity 


enters into the problem in a fluid which goes into motion 
from a state of rest anc on its effects toward preventing 

the existence of a velocity potential," Académie des Sciences 
Comptes Rendus, XC, (13880), 736-739; 967-953. 
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The inertial terms are given by 
YX ርነ =. E 0,0) 

ot 


The rate of deformation tensor is 


1 ow 
A. (14) = | o o e 


jio o 


and the exira stress tensor is 


Se i O O Due o o I! 
tk e (፦ቱ+?›)|] le ı ei + — lo 00 
à ۴ G O | /人 2 ( © © 


with 


Ky -  )16( ۰ 
The Characteristic Kelvin and Poynting effects of the non- 
Newtonian fluid theory are observed. If the additional 
normal stress be wanting in the vertical direction, the 
boundary wili tend to move up or down according to the 
sign of K,. The physical effect of the additional normal 
stress in the direction of fluid motion is not as easily: 


interpreted. We shall discuss this further below. 


Equations oi motion. - Recalling that ጁ (3) = (zt), the 


equatiens of motion take the form 
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(9-1) اس 
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O = Z-|-ptht+t ፳ (2 2271)‏ 


Solutions by Cases. 


سوریو خی وزو رت 
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Case I. ~ Classical case. Suppose Z = 元 “O, 7 = 2 At 
and £ is constant. Then (9-1) reduces to 


2 Ji 
= -ረዴጳ = - 2ው ተሥትሯ 
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Now from (9--2)o, ۳ is a function of X only. From (9-1); 


(9-2) 


e = O,.X% t Oa 


where Ou and A, are constants. Boussinesq^ Studied the case 


where Oz O and gave, as a solution to the flow equation, 


w= kt i- 5.1, ዕ- ዌ5) سو‎ (9-3) 


Lie 


Case IIA. - The relationship between assumptions con- 


wbere 


. cerning the pressure and existence of flows of non-Newtonian 

fluids, as discussed previously, is demonstrated in the prob- 
f 

lem of Boussinesq when 7; is considered constant, in which 


case (9-1) becones 
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where, fron (9.4) a 


0 14 (3) + Å (st). (9-5) 


Supposing, now, that we assume that is & constant. Then 
from the explicit relationship for ዮ› since 5) is a 
function of L and € only, it must be the case that 


oe A(t) and A, 5 ACC 
oF 
hence 


Ww = z Alt) ^A, (¢) መ. 


that is, the velocity is a linear function of the distance 


delt from the boundary and therefore the second derivative van- 


ishes. The right hand side of (9-4), therefore vanishes, and 
م = س‎ 
= 
t 
whence 


u = At + Alz). መና 
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À comparison of (9.8) with the above yields 


A, lt) = kt 
A(t) = comal. E K 


As(z) = Kz 
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so that the velocity is of the form 
utzx,ti=kt+Kz. (9-7) 


The boundary conditions, applied to this form, restrict the 
velocity field to the trivial case, (4 identically zero. 

We have therefore shown that, under &ssumption of 
constant pressure and ይ constant, there can be no Bousginesq 
flow of a non-Newtonian Huid. 

Viguier” discussed the problem of Boussinesq in the 
case where the velocity gradients are not negligible, em- 
ploying, in his equations of motion, viscosity forces s odd 
functions of the velocity gradients 8 sarlier introduced” 
by himself. 

Equation (9-7) comes about £8 £ direct consequence of 
easuming that fe and P gre constant. Since (42,9) *O , it must 
be the case, as noted from the conclusion above regarding the 
{mpossibility of the flow, that ሬ/,#|፡ ፊያ . Row the sssumption 
of adherence to the boundary gives w (0,6) s@ from which must be 
concluded that the velocity is identically zero. I have inter- 
preted this result as indicating that the sssumption of the 
constancy of pressure should not be msde and remove this 


condition in the cases studied below. 


EE س‎ 


3g. Viguier, "Some thoughts on & problem of M.J. 
Boussinesq," Académie royale de Be ue Bulletion Clssse 
des sciences, 1950-1, 71-76. 


fa. Viguier, "The equation of the boundary layer in 
the case of higher order velocity gradients," Académie des 
Sciences Comptes Rendus, CCXXIV, (1947), 713. 
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Viguier retains the assumption of constancy of presara 


and removes 959 condition 54. adherence, ንን ን "the condicion 
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Bowever one aky attempt to give partial justification." 

| Case IIB. 3 If we assume that the pressure is a 61 1“ 
of the flow direction and time only, similar considerations 
yield precisely the same results as in Case IIA above. 

Case IIIA. - We are inevitably lec in any study of the 
problem of Boussinesq, from the point of view of non-Newtonian 
fluids, to assumption of the dependence of the pressure diety i= 
bution on the distance from the boundary. We ‘consider A» A 

e constant and assume P و‎ 

Then from (9-1), we obtain, 

سا د مر i--k-‏ 


e az , (9-8) 
۱ ፦ ነ» 
(55)  * ቆዕ 
We may take, as a solution for the first equation, that given 


by Boussinesq 


+ حب‎ 
A = kt 7 اچ‎ (1- ን e A (5.5) 


ሠ. 


where, as before 
a Ah 
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Differentiating (9-9) with respect to Z yields 
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which can be rewritten, after integration by parts, 
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from which the pressure may be pm 


2 ጭ : 2 
peu a سنا‎ - 200 lech +A, (t) (9-1 


We consider the growth properties of the pressure. It 


is clear that the pressure is constant downstream. At the 


boundary, 


blot = - መረ r An). 


Returning to (9-10), 
-Jo (2,0) _ 
E. +0 
and the other term in the square bracket is easily seen to 
vanish as Í — O . Hence pee) approaches ፌ(ሮ) as one would 
expect. An additional initisl condition is necessary to establish 
the function / — uniquely. 


Now, as Z increases without bound 


wef e eto iso 
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and therefore plat) approaches ø (t) . 


in 


Finally, as Z increases, since [edn - SCH 
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p(z,t) جه‎ 9 


Case IIIB. - Supposing now that the pressure is a func- 


tion of X , Z , and É. Then 
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Ju M x ጋዶ 2 
eye - C^ 7 — Sx * "32 
(9-11) 

ነ — ጋሬ! 
and 如 = و‎ ٠٠ + (x,t). 

dÉ 

From (9-11),, since م‎ is independent of X 5x cannot be 


a function of X . Therefore, from (9-11)9, 


and therefore 
ቀ c, E) = x Alt) + A. (€). 


Hence we obtain, from (9-11), 


Yu 
سه‎ tAlt) = R+ 55 


with the integral e 
- a — 
2 Kê ۱ -4\u-#e "AN -ታ| ሐው 
we 
Now since {(0,€)20, and the term in the bracket vanishes 
at ZsO, it must be the case that fade vanishes identi- 
cally and therefore A, vanishes identically. 

We point out, finally that in the completely general 
case when fr = Ja) , from consideration of the general 
explicit form of the pressure, we see once more that the 
pressure must be a function of the distance from the boundary. 

In conclusion, we have demonstrated that the second 


order theory of the flow of an incompressible non-Newtonian 


fluid demands, in the problem of Boussinesq, that the 
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pressure distribution must be of Specific functional de- 
pendence. It cannot be considered constant without forcing 
the classical case or a trivial case and, further, the pres- 


Sure cannot be a function of the flow direction, 
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X. THE DRAINING OF A VERTICAL PLATE 


General. - We solve the equations of flow of a flat plate 
of infinite width which is drawn from an infinite fluid and 
compute the thickness of the fluid film on the surfaces of 
the plate. We assume a condition of steady state has been 
reached, that the fluid is incompressible and that the vel- 
ocity is vertical and independent of the height. 

This problem was first considered by Goucher and 
Ward? and an equivalent problem studied by Jeffreys”. 

We consider a rectangular cartesian coordinate system 
with the Y-axis horizontal and perpendicular to the vert- 
ical plate? the Z, -axis is vertical and parallel to the 
plate, so that when EM , the thickness of the film, uf: 0 
and "یا‎ > 0 and that ya V, a constant (i.e. , the vertical 
velocity of the plate), when 4*0 , UF being a function of 


d only. We disregard the curvature of the surface in the 


47r plane and assume the acceleration of the fluid is 
negligible. 
Undor the above assumptions, the velocity vector is 


given by, 





l.S. Goucher and H. Ward," A problem in viscosity: 
the thickness of liquid films on solid surfaces under dynamic 
conditions," Philosophical Magazine and Journal of Science, 
XLIV, (1922), = - =. 


Sa Jeffreys, "The draining of a vertical plate," 


Proceedings of the Cambridge Philosophical Society, XXVI, 


This is equivalent to assuming that there are no 
normal stresses on the free surface. 


—— 
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X (i) = (0, 0, ww). 
The boundary corrditions are: 


aur (T) s 
w'(T) *» 0 
ar (0) = Y 


the external foiwc$ fisld is 
FG) 2 Coa =g) 


and the rate 0:1 deformation tensor is 


4) 0 ٥ 
ዲ(፥ከ) ። lo o qu 


3 ትሠ' ጋ 


The stress tens: takes thie form 
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where 
Tas fl oo ` I=-zw” A “/ወሠ. 


and the Poynting effect is noted. An additional stress along 


the ክነ is soted; if this stress be wanting, the thickness 


of the film wil] tend. to increase. 


, Equations of motion. - The equations of motion take the form 
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E 
o = sf 
O = 51 - لد + د‎ ur”) (19-15, 


gd‏ د اس یون وم 

From (10-1),, if 加 is constant, then Fa ሾሆ" is constant. 

But Fa wi? is a power series in wr, and therefore it must 
be the case that wi ís constant. Since ur’ vanishes when 
st s ig’ must be identically zero, hence WM must be constanti. 
Since YW= 0 when 4*7 , uf must also vanish identically and thx 
flow is trivial. We have demonstrated that, in the draining 
of a plate being removed from o non-Newtonian fluid, tħe preas- 


Sure cannot be assumed constant, as in the classical case. 


The solution. - Since the pressure is independent of X* we 


have from (10-1)9 

和 三 5 w' *AG).- (16-8) 
Integrating (10-1)3 with respect to Z , we obtain 

P=2 E (k ur’) -ረዳ + A.(y). | (10-55 
Comparing (10-2) and (10-3), we may write 

P = (^i- 63) Z A ها و لل‎ + cost. 


and A 


ho) 2 ۸, 2 روم‎ 


the choice of A, following from the absence of normal stresses. 
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By integration and, from the boundary condition for 


LJ! , we obtain 
> f, ur’ = A, (4 -T). (10-4) 


As in chapter V, we revert (10-4), obtaining 


20, A8?A; — orano aaie 
w= 251957) ፥ ማመ 3! ፥ 


which, upon integrating and substituting (3-5)o, 


— ABA (st EES 
LU = په‎ 1 + TP dh (10.5) 
When =O 
8 AT? a. ABA 4 
uy = y= E e XI md ት . “ 95 5 ቀ . © ® (10.6) 


We may therefore combine )10-5( and (10-6) and write 
5 A CH ne 


3 AJ 
ተ ሚዎ (4-44 T*64 T - 44 T) — 


The first approximation, 


«a sans Geer 
is the solution of Goucher and Ward. 

We now seek numerical values for the thickness, 7 , 
in the second approximation, in order to compare these values 
with those reported experimentally by Goucher and Ward. To 
this end, we take only the first two terms on the right in 


the expansion given by (10-6), and may therefore write 
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$ A * As - 0۸ 
a." T مره‎ Y = 0 (10-7) 


where © = a AS and is & constant dependent only on the 


fluid much the same 85 JA . Me see that 4,20 implies that 
-0 and 


لا و A‏ 


we write the expression for the square of the thickness 


3 842 BA 
A ባፈ A € (10-8) 


In order to compare the experimentally measured thicknesses 


with those given in (10-8), we take 
pm=.25, A =/420 


in c.g.s. units, which fall within the range of the measure- 
ments of Goucher and Ward. 
Using & point from the experimental results, we then 
MT 


compute an value for in c.g.S. units (dimension T )m 


- 3 
中 is negative and in the interval 
-10< €«10% ¢-] 


with a value of -8 6 ም yielding a curve which compares 
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very well with observation. We rewrite «¢ for ዌ in (10-8) 


7 تا 5 د01 


From knowledge of the dependence of the thickness, | , on the 
velocity, Y » in the classical case, and from the available 
experimental results, we know that the thickness increases 
with an incrense in the velocity. Therefore the expression 
in (10-9) with the negative sign is that which yields the 
solution to the physical problen. 

The following table shows comparative values in the 


classical, non-linear and experimental cases: 


VELOCITY THICKNESS 
Linear non-linear — experimental 

0,00 . DADO . 0000 . 0000 
0.16 .0075 .0076 .0078 
0.64 .0150 ‚0151 «0165 
1.44 0225 0229 .0255 
2.56 . 0300 .0311 .0345 
4.00 . 75 .0398 , 0438 
5.76 ‚0450 . 0494 ` .0546 
7,84 ¿0525 . 0610 .0669 
10.24 .0600 ‚0820 .0822 


The above values of T are plotted as 2 function of ۷ in 

Piles El ` 
Now, in the classical case, the thickness increases 

| without bound, ag the velocity increases. Tatu tive, e 

would expect that, for a given fluid, i,e,, for a given 

$, م‎ + and 4 › there probably is a velocity beyond which 

the assumptions regarding laminar flow and boundary condi- 


‘tions at y=T are no longer valid. This is evident in the 


non-linear theory. 
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Tie expression in (10-9) no longer has physical meaning 


Vt the point 


x^ tthe cieht 


2 
SA, 
where the velocity is equal to the expression 


in (10-10), we obtain 


the maximum value that can be taken on by Y, so 


taxt the maximum value attainable for the thickness is given 


ty for a particular fluid 
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APPENDIX 


DERIVATION OF THE ZQUATIONS OF HYDRODYNAMICS? 


1. Equation of continuity. - The equation of continuity ex- 
presses the conservation of mass. Suppose that a sufficiently 
Smooth surface, S , includes a volume of fluid mass, y. 
We define the mass, M, 

ነሳ 3 ነር” 
vhere j indicates integration over the volume V ; £ is the 
fluid density and AY an element of volume. We wish to ex- 
press the principle that mass is neither created nor destroye:. 
Hence the time rate of flow of mass from a given stationary 
volume must vanisb, except for any mass which is added or sub. 
tracted from the volume. 


4 
We let - be the rate of increase of mass. The 


N 


local mass flux acrose an element of area (directed outwardly) 
As; per unit time is given by ex As; where X’ is the velc- 


city vector. Therefore to express conservation of mass, we 


ኣ 


write, 
dm . -$ x* d$; 
dt 5 e 
whore $ jndicates integration over the surface 9. Hence 


dr len ኑን መን = ር2 


Assuming (Y^ and f continuously differentiable, ar. 


f 
applying Green's Theorem, we obtain 


\ ES + ,(۵ع)‎ 2V=0 
۷ 


1 Following the manner of Truesdell in lectures on 
Hydrodynamice at Indiana University, lst semester, 1951-19582. 


TE 
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whence we obtain the Eulerian Equation of Continuity, 
ee NN. pip - = O 
A ፡ድዶ”/ዞቆ 
2. Cauchy's Laws of Motion (1821, 1828) From an axiomatic 
approach to continuum mechanics, we have that the force and 
moment acting on a body consists of three parts: 

1. <A force and moment equal to resultant force aud 


ጄ 
moment, respectively, of a certain vector field, Z'(,) , de- 


fined on the boundary, 3 , of a body; thus 
j ዖዳ 
= ab E. ቋ5 Li = Erjk XC) AS 
3 


where رود ع‎ / if ¿jk are an even permutation of /,2,3 
ca: ~/ if an odd permutation or @ if z,j,A are not all dif- 


forent. We call — the stress vector and shall see that 
M_ 


dimension Ets, is FEA 


11, A force and moment equal to the resultant force 
anc momentum respectiivery of a certain density, 


f \ ~ ጥኑ › x; e > t) 


where P , Q are two points. Then 
' rk 0 
lk dam e = e 


where ٢ سا‎ dung) over the whole space represents 
tihe effect of Sea on a point P as Q passes through all 
“paints of the space. 


111. A force and moment equal to the resultant force 
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and moment respectively of a certain moment density T CA xt, +) 


plus a certain noment density Åe Ott) ; then 
i , SA 
F=f fh dew d = [ei e the.) daw 


and fhe fr rk is the extraneous force field. 


We therefore have the axioms of continuum mechanics 
p-$ Ej 85 + A m x ALL 
$ k 1 k 
> T * | g Pd 
P & Esje Em, AS d ior Xf fi) he = J 
* 


O 
t 1 
where ai and Ber are the linear and angular momentum 


respectively. The subseript (») refers to the side of a 


፥ 
particular surface, 9 , on which ያሬ acts. 


3. Cauchy's Laws of Motion. 
i. Cauchy's First Theorem Let $ be any surface 
(assumed regular enough to allow application of Stokes’ aind 
Green's Theorems), bon) and Ein) be the vector densitime 
| on the + side and -- side of the surface, § , each contiin- 
uous on its own side; then if po x* is of bounded variation, 
too 5 t 
Proof: This theorem is proved from the equatiom 
for p applied to a surface on each sideof S . 
Taking as a bound on the integral of ን፦ኑሯ , ind 


allowing the surfaces to shrink to S , we obtain 
3 pt 
A Ba * | Line As = 0 
5” st 


whence from the assumed continuity, 


? 2 
tay * t (n) 
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ii. Cauchy's Second Theorem Using the expression 


y ጌ ። 
for P we show the existence of a stress tensor, defined by 


` Lt 
Tin 5 * A 
where A; is the unit normal to 2 surface (directed outward). 


The use of the stress tensor allows us to describe the stresses 


on a surface without having to speak of the stress direction. 


ili.  Cauchy's Laws of Motion. - We consider only the 


continuously differentiable case. From the expression for Pp 


$ tints +8 Im = O. 
3 v 


Now, 


Ó tin as 3 ቁ en; 85 = h “ds; = | t۷ 
5 5 3 Y 


by Green's Theorem. Upon substituting ዲላላ ። eav, we obtain 


\ je; rel SI 8ኛ = 0 
. ۱ | 
of t'ieg(f'-xX) 20 
Cauchy's First Law of Motion. 
From the expression for SÉ Wwe can similarily show 
that ek &* , of; = O , Cauchy's Second Law of Motion. 
Now if E ; the extraneous moment of momentum 


vanishes, then Ltd war ‚ which we henceforth assume. 


4. The Energy Equatuion. - We define the kinetic energy, K., 


ረ= £ fex xi 
V 


i 

" "UI 
EE 
ur — 


ed 


SAS 





E 
We further consider the internal energy ራድ of a body, and 
assume that it is a continuous function of mass. We may 
then consider an integral form Ez sim where € is the in- 
ternal energy density. E is changed by two specific pro- 
cesses, (1) mechanical work and (2) one TG causes 


(e.g. radiation, heat conduction, etc.) 


Principle of Energy. - Let there be a heat flux, 


j » (analogue of the stress vector) such that we may write 
4 e E £s V pr y Bi: سه‎ > Ds. 
KE = | 61 ደሃ ቶፉ Ki Rij -P F° AS 
Vu 3 ሠ 


that is, the total time rate of change of energy is equal 

to the mechanical work performed plus the work done due to 
the stress tensor. p gias: is so defined as to balance the 
equation. We note that the interconvertability of heat and 
work is already determined by this equation because of the 


fact that st is written in mechanical units, i.e. 


dimension & کم ۹625۷ د‎ .M. 
time area Eé 


Assuming all functions involved are continuously differenti- 
able, substituting the expression for A and € and applying 


the divergence theorem, we obtain 
EA 5 4 ቃ yb y on 
\ | x (e i ږې په‎ | aV AU - "Xj +3° 
9 k - 


Applying Cauchy's first law, and from continuity assumptions, 


we have 
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Defining the rate of deformation tensor as usual, 


ቓ1።፡ ጆን)‏ = ږل 


and recaliing the symmetry of the stress tensor, we have 
PEt - Ba the Fourier-Kirchoff-C. Neumann equation. 
Thus £94; ; y the stress power, is the rate at which stress 
does work per unit volume. 

Assuming as usual for a homogeneous fluid that £ is given 
by & = €»), where V ^47 and 7 ís a local statistical 
parameter (entropy), we define 77 = “#«923 , the pressure 
and temperature respectively; then ARE 2=>27” ያሇ 四 ያት. 

We now define the extra stress tensor vi by 
vi = #5 + E , where ያ is any scalar. After some al- 


gebraic manipulation, wo obtain 

ፇይ =- 7 ፈ + گم هن‎ - f 
Hence 

PON = priv f ZH # ሆ* ከዕ SCH 
From the equation of continuity 

div xk = ዲኒ = ገሆ 
whence 

007 = አረን jy + (T-p) Y 
For an بر مه‎ Mg fluid Y 20, and 


por = vide - 25: , 
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